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ABSTRACT 

We present the largest set of N-body and hydrodynamical simulations to date for cos- 
mological models featuring a direct interaction between the Dark Energy (DE) scalar field, 
responsible of the observed cosmic acceleration, and the Cold Dark Matter (CDM) fluid. With 
respect to previous works, our simulations considerably extend the statistical significance of 
the simulated volume and cover a wider range of different realizations of the interacting DE 
scenario, including the recently proposed bouncing coupled DE model. Furthermore, all the 
simulations are normalized in order to be consistent with the present bounds on the amplitude 
of density perturbations at last scattering, thereby providing the first realistic determination 
of the effects of a DE coupling for cosmological growth histories fully compatible with the 
latest Cosmic Microwave Background data. As a first basic analysis, we have studied the 
impact of the coupling on the nonlinear matter power spectrum and on the bias between 
the CDM and baryon distributions, as a function of redshift and scale. For the former, we 
have addressed the issue of the degeneracy between the effects of the coupling and other 
standard cosmological parameters, as e.g. as, showing how the redshift evolution of the 
linear amplitude or the scale dependence of the nonlinear power spectrum might provide 
a way to break the degeneracy. For the latter, instead, we have computed the redshift and 
scale dependence of the bias in all our different models showing how a growing coupling 
or a bouncing coupled DE scenario provide much stronger effects with respect to constant 
coupling models. Furthermore, we discuss the main features imprinted by the DE interactions 
on the Halo and Subhalo Mass Functions. We refer to this vast numerical initiative as 
the COupled Dark Energy Cosmological Simulations project, or CoDECS, and we hereby 
release all the CoDECS outputs for public use through a dedicated web database, providing 
information on how to access and interpret the data. 
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1 INTRODUCTION 



The observational evidence in favor of a relatively recent acceler- 
ation of the expansion of the Universe has been progressively re- 
inforced by a large number of independent probes (see e.g. 



et al. 



2005 



Astier 



20061 Kowalski et al. 2008 , Per cival et al.|2001| |Cole et al 



|Percival et al.|2010||Reid"et al.|2010) over more than a decade 



since the first discovery of the dimming of distant Type la Super- 
novae (Riess et al. 1998, Perlmutter et al. 1999). The fundamental 
origin of such acceleration - generically referred to as the "Dark 
Energy" - represents one of the central open problems in modern 
cosmology, and is the main target of a wide range of ambitious 
observational initiatives planned for the upcoming future such as 
PanStarrs l |Kaiser et al.|[2002) , HETDEX l |Hill et al.||2008l >, DES 



I Abbott et al.|2005f , LSST |lvezic et al.|2008} and Euclit^] |Lau 
|reijs et al.|20TT) . 

The standard ACDM cosmological model accounts for the ob- 
served Dark Energy by means of a cosmological constant A, i.e. a 
completely homogeneous form of energy whose density does not 
dilute with the expansion of the universe. Such a simple prescrip- 
tion is sufficient to fit the vast majority of presently available data, 
but requires an unnatural level of fine-tuning of the cosmological 
constant value, giving rise to the so-called "fine-tuning" and "co- 
incidence" problems of the standard ACDM cosmology (see eg. 
|Weinberg|1989"l >. For this reason, alternative Dark Energy scenarios 
have been proposed, including models based on dynamical scalar 

1 http://www.euclid-ec.org 
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fields jWetterich| 1988j |Ratra & Peebles! 1988| |Armendariz-Picon| 
|et al.|200T] > or on modifications of General Relativity at cosmolog- 
ical scales (see e.g. |Hu & Sawicki|2007t . 

Furthermore, a wide range of astrophysical data at intermedi- 
ate and small scales seem to show possible deviations from the ex- 
pectations of the ACDM model, thereby providing additional moti- 
vation for the exploration of alternative scenarios. Such deviations 
range from the lack of luminous satellites in CDM halos (the so- 
called "satellite problem", Na varro et al.||199"6), to the ob served 
low baryon fraction in galaxy clu sters (see e.g. | Ettori|2003]|Allen 



et al.|2004||Vikhiinin et al.|2006||LaRoque et al.|2006[|McCarthy 



et al.|2007 1, to the shallow observed CDM density profiles in dwarf 



galaxies (the so-called "cusp-core" problem, |Moore 1994 



Flores 



& Primack[ l994 > Sim on et al.|2003} , in spiral galaxies ( Navarro & 
Steinmetz|2000||salucci & Burkert|2000||Salucci|200l||Binney &| 
Evans|2001), and in galaxy clusters ( |S and et al.|2002||2004||New-| 

man et al. 2009), to the anomalous detection of massive clusters 



at hig h redshift <|Mullis et al.|2005| |Bremer et al.|2006[|Jee et al 



2009[|Rosati et ai.|2009||Brodwin et al.|2010||Jee et al.|201 l[|Foley 
et al.|201l| , and to the high velocities detected in the large-scale 
bulk motion of galaxies (e.g. Watkins et al. 2009J or in systems of 
colliding galaxy clusters as the "Bullet Cluster" jLee & Komatsu| 
[20T0l|Lee|20T0l >. 

Most of these anomalous observations involve highly nonlin- 
ear environments where complex astrophysical processes play a 
major role in shaping the evolution of structures and of their di- 
rect observable properties. It is therefore reasonable to investigate 
a possible astrophysical origin of such deviations from the expec- 
tations of an idealized ACDM cosmology by refining our modeling 
of all the relevant astrophysical processes in place at the scales of 
interest for these different observations, as e.g. gas cooling, star for- 
mation, and feedback mechanisms from supernovae explosions in 
galaxies or from Active Galactic Nuclei in galaxy clusters. Never- 
theless, such deviations could also represent the signature of the un- 
derlying cosmological model, thereby providing a unique handle to 
investigate alternative scenarios to the standard ACDM cosmology 
and to constrain the properties of the dark sector of the universe. 

In both cases, however, due to the high nonlinearity of the 
systems under investigation, it is necessary to rely on complex and 
computationally demanding numerical simulations where the spe- 
cific astrophysical or cosmological models to be studied have to 
be properly and efficiently implemented in order to provide reli- 
able predictions to be compared with real astrophysical data. In 
this respect, if on one side a great effort has been made in order 
to refine the theoretical modeling and the numerical implementa- 
tion of complex astrophysical processes in cosmological N-body 
codes (see e.g. |Springel & Hernquist|2003||Maio et al.|2007| \st\ 
jacki et al. 2007 ; Springel 201 1 1, not as much deal of work has been 
put into the modification of standard N-body algorithms in order to 
account for the nonlinear evolution of alternative cosmological sce- 
narios. In the last years, some first steps have been moved in this 
direction by a few numerical studies of non-standard cosmological 
scenarios as e.g. Modified Gravity models (in the form of f(R) the- 



ories, [ Oyaizu 2008; Oyaizu et al. 2008l |Schmidt et al.|2009l|Zhao| 



et al.|201 1), Warm Dark Matter models (see e.g. |Zavala et al. 2009| 
Lovell et al. 2012; Viel et al. 2011), Extended Quintessence mod- 
els jDe Boni et al.|20TTj|Li et al.|2011) , or Coupled Dark Energy 
models ( |Macci6 et al.|2004||Baldi et al.|2010|p & Barrow|20"TT] >. 

In particular, the latter have been shown to determine signif- 
icant deviations from the expectations of ACDM in e.g. the pre- 
dicted baryon fraction and concentration of galaxy and cluster halos 
( |Baldi et al.|2010| |Baldi|201l"b} , the predicted transmitted flux of 



Lyman-Q systems (Baldi & Viel 2010), the number counts of mas- 
sive galaxy clusters as a function of redshift ( Mainini & Bonometto 
|2006||Baldi & Pettorino(20TTl|Baldi|2012||Tarrant et al.|2012) , and 
the distribution of galaxies within a galaxy cluster main halo (Baldi] 
|etal.|2011| ). 

Coupled Dark Energy models have been first proposed by 
|Wetterich|(T995l l as a possible way to alleviate the cosmological 
"fine-tuning" problem, and subsequently studied by several authors 
(see e.g. |Amendola|2000||2004||Farrar & Peebles|2004{|Mainini & 
Bonometto|2006| Farrar & Rosen 2007 |Pettorino &~B accigalupi 



2008; Amendola et al. 2008; Koyama et al. 2009, Caldera-Cabral 
eraTi|2009; Honorez et al.|12010 Baldi||201 lbl 12012] 



These are 

based on the dynamical evolution of a classical scalar field cj> that 



plays the role of the Dark Energy, and feature a direct interac- 
tion between this Dark Energy scalar field and the Dark Matter 
fluid, that exchange energy-momentum during the cosmic evolu- 
tion. Such interaction gives rise to a "fifth-force" between Dark 
Matter particles - mediated by the Dark Energy scalar field - 
that significantly modifies the gravitational instability processes 
through which cosmic structures develop, both in the linear and 
in the nonlinear regimes. A detailed investigation of the nonlinear 
evolution of this particular class of alternative cosmological mod- 
els is therefore necessary to provide a quantitative determination 
of the deviations expected from the standard ACDM predictions 
in several observational contexts, and to assess the viability of the 
Coupled Dark Energy scenario as a possible alternative to the stan- 
dard model. 

In this work we will present the largest and most accurate N- 
body simulations to date for Coupled Dark Energy cosmologies 
in terms of simulated volume, numerical resolution, and range of 
models covered in the numerical sample. These include both colli- 
sionless runs at large scales and adiabatic hydrodynamical simula- 
tions at small scales for five different Coupled Dark Energy scenar- 
ios, besides the standard fiducial ACDM cosmology. The various 
Coupled Dark Energy models include constant coupling models 
(Amendola 2000), variable coupling models (Baldi 201 fij}, and the 
recently proposed Bouncing Coupled dark Energy scenario jBaldi| 
|2012[ >. All the models share the same set of cosmological parame- 
ters at the present time, and the same amplitude of density perturba- 
tions at the redshift of the last scattering surface (zqmb ~ 1100), 
both consistent with the latest results from the WMAP satellite ( Ko- 
|matsu et al.|201 1[ ). The present set of cosmological simulations is 
therefore fully compatible with CMB constraints and provides a 
detailed picture of the evolution of the universe from the last scat- 
tering surface to the present epoch within five different Coupled 
Dark Energy scenarios. 

We will therefore illustrate the details of this vast numerical 
study that goes under the name of the "CoDECS project", where the 
acronym CoDECS stands for (COupled Dark Energy Cosmological 
Simulations). All the simulations outputs - that include, besides 
the raw snapshot files, also halo and subhalo catalogs and the 
nonlinear matter power spectra - are made publicly available 
through a dedicated web database together with the background 
and linear perturbations evolution computed for each different 
model. 

The paper is organized as follows. In Section|2]we review the 
basic equations and main features of Couped Dark Energy mod- 
els both concerning their background and liner perturbations evo- 
lution. In Section [3] we introduce our set of N-body simulations, 
describing their numerical features and the procedures adopted to 
generate initial conditions, and comparing the visual outcomes of 
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a few sample simulations for different Coupled Dark Energy sce- 
narios. In Section[4]we present some basic results derived from the 
analysis of our large-scale collisionless simulations, focusing on 
the time evolution of the nonlinear matter power spectrum and on 
the offset between the baryonic and the Cold Dark Matter overden- 
sities. We also provide a few examples of the possible applications 
of the CoDECS halo catalogs by computing the Halo Mass Function 
and the Subhalo Mass Function for the different cosmologies under 
investigation, and highlighting the main characteristic features im- 
printed by the Dark Eenergy interaction on these observable quan- 
tities. Finally, in Section|5]we draw our conlusions. 

Furthermore, in the Appendix we provide the details of the 
public database where the simulations outputs are stored and acces- 
sible for public use, briefly describing the contents of the database 
and how different information are organized in different files. High- 
resolution versions of the images displayed in this paper can also 
be downloaded from the same website. 



2 INTERACTING DARK ENERGY MODELS 

The CoDECS project is aimed at providing publicly available data 
from large N-body simulations for a significant number of interact- 
ing Dark Energy (DE) cosmological models. Interacting DE sce- 
narios have been extensively studied in the last decade as a possi- 
ble alternative to the standard ACDM cosmology as they feature a 
dynamical evolution of the DE density that provides a way to par- 
tially address the fine-tuning problems that affect the cosmological 
constant (see e.g. [Weinberg 1989} and to simultaneously account 
for an accelerated expansion of the universe at late times. Starting 
from the initial proposal of Wetterich ( 1995), this class of cosmo- 
logical models has been investigated by several authors ( Amendola 
2000) [2004) |Pettorino & Baccig aTuppOQ8[| Amendola et al.|2008| 
Caldera-Cabral et al.|2 009, Koya ma et al.|2009j > to which we refer 
for a detailed derivation of the main equations and for an extensive 
discussion of the main features of coupled DE (cDE, hereafter) cos- 
mologies. 

In the present Section we will only briefly review the basic 
equations that characterize cDE scenarios and we will present the 
specific models that are included in our suite of N-body simula- 
tions, providing a direct comparison of their background and linear 
perturbations evolution. 



2.1 Background evolution 

We consider a set of flat cosmological models that include CDM, 
baryons, radiation and a classical DE scalar field (f> (subscripts 
c, b, r, and <f>, respectively). For any given set of cosmological 
parameters Q, r , Q,b, Q<t>, flatness is ensured by imposing Q c — 
1 — fl r — fij, — In particular, for all the models included in the 
CoDECS suite we consider cosmological parameters in agreement 
with the "WMAP7 only Maximum Likelihood" results of |Komatsu] 
[etaLlpOTT) , that are listed in Table[T] 

The background dynamic equations for the different cosmic 



Parameter 


Value 


H 


70.3 km s" 1 Mpc" 1 


^CDM 


0.226 


^DE 


0.729 


A s 


2.42 X 10~ 9 


n b 


0.0451 




0.966 



Table 1. The set of cosmological parameters at z = assumed for all the 
models included in the CoDECS project, consistent with the latest results 
of the WMAP collaboration for CMB data alone iKom atsu et al.|201 1) . 



components are given by 

4> + 3H<j> + = 



p c + 3Hp c 
Pb + 3Hpb 

p r + 4:Hp r 



5«< 



3ir = 



(Pr + Pc + Pb + P4>) 



(1) 

(2) 

(3) 
(4) 

(5) 



where the source terms at the right hand side of Eqs. \ \\-\ encode 
the interaction between the DE scalar field cj> and CDM particles. In 
Eqs. \l\2\ , and throughout the paper, the scalar field <j> is expressed 
in units of the reduced Planck mass Mpi = 1 /v&kG and an over- 
dot denotes a derivative with respect to the cosmic time t. The cou- 
pling function /3 C {4>) sets the strength of the interaction while the 
sign of the quantity (j>/3 c (<j)) determines the direction of the energy- 
momentum flow betw een the two components. With the convention 
assumed in Eqs. |l|2| a positive combination <^/3 c (<j!>) > corre- 
sponds to a transfer of energy-momentum from CDM to DE, while 
the opposite trend occurs for negative values of <j>/3 c (<f>). Such en- 
ergy transfer implies a time variation of the CDM particle mass, 
according to the equation: 



m c (z) 



m c ,o ■ e 



- ! Po(4>)d<t> 



(6) 



which can be obtained by integrating Eq. j2|, that implies rh c < 
for (j>/3 c ((f>) > and rh c > for 4>fi c (<f>) < 0. As shown by 
Eqs. ( |3|4| >, baryons and radiation are always uncoupled from the 
DE scalar field, such that pb oc a~ 3 and p r oc a~ 4 for all the 
cosmologies considered in the present work, and Eq. (|5| closes the 
system by imposing flatness. 

In the range of models included in the CoDECS project we will 
consider two possible choices for the scalar field self-interaction 
potential V ((/>), namely an exponential potential (Wetterich 1988): 



V(<t>) = Ae~ a 
and a SUGRA potential ( |Brax & Martin|1999| >: 

V{<f>) = M' 



-a S' s /2 



(7) 



(8) 



The former is known to provide stable scaling solutions for the 
scalar field that are almost independent from the field's initial con- 
ditions even for an uncoupled system (i.e. for Pc{<j>) = 0, see e.g. 
Ferreira & Joyce|1998| >. For the case of coupled DE (Pc{4>) / 0), 
instead, it features a transient early DE scaling solution (known as 
the (^-Matter Dominated Epoch, or <^>MDE) followed by a late-time 
accelerated attractor (see e.g. Amendola 2000, 2004, Baldi 201 lb I. 
The latter potential form, instead, naturally arises in the context of 
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Model 


Potential 




a 


Po 


01 


Scalar field 
normalization 


Potential 
normalization 


w^(z = 0) 


A s (zcmb) 


a$(z = 0) 


ACDM 


V(<f>) = A 












A = 0.0219 


-1.0 


2.42 x 10~ 9 


0.809 


EXP001 


V(4>) = Ae~ 




0.08 


0.05 





<f>(z = 0) = 


A = 0.0218 


-0.997 


2.42 x 10~ 9 


0.825 


EXP002 


V{4>) = Ae~ 


■Ct(j) 


0.08 


0.1 





<j>(z = 0) = 


A = 0.0218 


-0.995 


2.42 x 10~ 9 


0.875 


EXP003 


V(4>) = Ae~ 




0.08 


0.15 





<f>(z = 0) = 


A = 0.0218 


-0.992 


2.42 X 10 -9 


0.967 


EXP008e3 


V(<j>) = Ae~ 




0.08 


0.4 


3 


<j)(z = 0) = 


A = 0.0217 


-0.982 


2.42 X 10- 9 


0.895 


SUGRA003 


V(<f>) = A4>~ 


-a £ <l> 2 /2 


2.15 


-0.15 





4>(z — > oo) = y/a 


A = 0.0202 


-0.901 


2.42 X 10- 9 


0.806 



Table 2. The list of cosmological models considered in the CoDECS project and their specific parameters. The scalar field is normalized to be zero at the 
present time for all the models except the bouncing cDE scenario, for which the normalization is set in the very early universe by placing the field at rest in its 
potential minimum (see Baldi 2012 for further details). All the models have the same amplitude of scalar perturbations at zcmb ~ 1100, as shown by the 
common value of the amplitude As, but have very different values of erg at z = 0, again with the sole exception of the bouncing cDE model SUGRA003. 



supersymmetric theories of gravity and provides a viable cosmo- 
logical evolution with a late-time accelerated attractor even in the 
absence of any coupling, while for the case of a cDE scenario it 
has been recently shown ( jTarrant et al. 2012 Baldi 2012 1 to have a 
significant impact on the evolution of density perturbations. In par- 
ticular, Baldi (2012) has extensively investigated the same SUGRA 
cDE model presented in this work, showing that it determines a 
"bounce" of the DE equation of state ui^ on the cosmological con- 
stant barrier = — 1 at relatively recent times, which has a signif- 
icant impact on the expected number density of massive CDM ha- 
los at different cosmic epochs. Due to this peculiar feature this type 
of models has been dubbed the "bouncing" cDE scenario ( |Baldi| 
|20T2l >. 

While in the case of the exponential potential models the 
scalar field always rolls down the potential with positive velocity 
4> > 0, reaching the normalization value <f> = at the present 
time, the bouncing cDE scenario based on the SUGRA potential 
is assumed to start with the scalar field at rest in the potential 
minimum in the early universe, and the dynamics induced by the 
coupling allows for an inversion of the direction of motion of the 
field during the cosmic expansion history. The peculiar dynamics 
of the bouncing cDE scenario has been described in full detail 
by Baldi (20121, to which we refer for a more comprehensive 
treatment of the basics of this model. 

To complete the definition of the range of models considered 
in the CoDECS project we need to specify the coupling function 
Pc(4>), for which we will assume the exponential form proposed by 
|Amendola|p004t ; [Baldi1j201 lb| >: 

= p e M (9) 

and we will consider in our analysis both the case of a constant 
coupling (/3i = 0) and of an exponentially growing coupling (Pi > 
0). 

All the cosmological models included in the CoDECS suite, 
with the relative parameters and normalizations, are summarized in 
Tabled 

For each of our cosmologies we compute the background evo- 
lution by numerically solving the system of coupled differential 
equations \\\5) . For the ACDM and the exponential potential mod- 
els EXP001-EXP008e3 we integrate the equations backwards in 
time starting from the desired set of cosmological parameters at 
z = 0, by applying the integration procedure described by |Baldi| 
( |2011b ). For the SUGRA cDE model, instead, since the normaliza- 
tion of the field is defined at high redshifts, we integrate the system 
forward in time with a trial-and-error procedure until the desired 
set of cosmological parameters at z = is obtained. 

The evolution of the DE equation of state parameter wj,, of 



the Hubble function H(z), of the CDM mass m c (z)/m c ,o, and of 
the DE fractional density Q$ are shown in Fig. [T| in panels A, B, 
C, and D, respectively. 

As mentioned above, the inversion of the field motion in the 
SUGRA cDE model is reflected in a bounce of the equation of state 
w,j, on the so called cosmological constant barrier (w^ = —1) at 
Zinv = 6.8, as clearly shown in Fig.[T|A. As one can see from the 
plot, all the cDE models have a stiff equation of state at high red- 
shifts while they all end up at z — with values of w$ that are 
very close to —1 (see also Table|2j, with the only exception of the 
bouncing cDE model that due to the bounce has a somewhat larger 
equation of state parameter at the present time (w$ = —0.91). As 
a consequence of its non trivial dynamics at low redshifts, the spe- 
cific bouncing cDE model considered here might therefore be chal- 
lenged by present observational bounds on the DE equation of state 
from low-redshift probes, as e.g. Supernovae la. As already clearly 
stated by |Baldi|pfJT2) , the SUGRA003 model should be taken as a 
phenomenological toy example of the effects that a bounce of the 
DE scalar field can imprint on structure formation processes, while 
a full parameter fitting (that goes beyond the scope of this work) 
would be required in order to select the best set of parameters for a 
generic bouncing cDE cosmology. 

The effects of the DE bounce on the background evolution of 
the universe are clearly visible also in Figs.[T]B and[T]C, where the 
Hubble function and the CDM mass evolution are plotted, respec- 
tively, as a function of redshift. The former plot shows a clear in- 
version of the deviation from the ACDM Hubble function in corre- 
spondence of the bounce for the SUGRA003 model; the magnitude 
of the deviation is of the order of ~ 6% at low redshifts and then 
slowly decreases towards very high redshifts. Such feature does not 
appear, as expected, in all the other cDE models based on the ex- 
ponential potential of Eq. j7). In the latter plot, instead, the effect 
of the bounce is shown by the inversion of the mass evolution for 
the SUGRA003 model, while all the other scenarios show a mono- 
tonic decrease in time of the CDM particle mass. It is important to 
notice here how the CDM mass in the SUGRA model has roughly 
the same value at zcmb ~ 1100 and at the present time. 

Finally, Fig.[T|D shows the evolution of the DE fractional den- 
sity as a function of the e-folding time In a for all the models 
considered in this work. The early DE nature of cDE models is 
clearly visible in this plot, where all the different cDE cosmologies 
show a non negligible DE fraction at high redshifts, which in any 
case never exceeds the percent level. The effect of the bounce in 
the SUGRA003 model is clearly visible also in this plot, with an 
evolution of 0,$ that significantly deviates both from ACDM and 
from all the other exponential cDE models at low redshifts. 
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Figure 1. Panel A: The equation of state parameter of DE as a function of redshift for all the models of the CoDECS project. The SUGRA003 model shows 
a "bounce" of the equation of state on the cosmological constant barrier = — 1 at z- mw ~ 6.8 while all the other models asymptotically tend to wj, = — 1 
at z — > 0. Panel B: The upper plot shows the Hubble function in km/s/Mpc for all the models under investigation, while the lower plot shows the ratio of the 
Hubble function to the standard ACDM case. The maximum deviation is of about 6% at z ~ 5 for the bouncing cDE model SUGRA003. Panel C: The CDM 
particle mass evolution as a function of redshift. All the models feature a monotonically decreasing mass, except the bouncing cDE model SUGRA003 that 
shows an inversion of the mass evolution in correspondence to the DE bounce. Panel D: The redshift evolution of the DE density parameter Q^. The constant 
coupling models show a constant plateau at high z (the so called </>-MDE regime), while the growing coupling model EXP()08e3 shows a slowly evolving DE 
fraction (the so called Growing- </>MDE regime, see |Baldi|201 lb} . The bouncing cDE model SUGRA003 behaves like a standard constant coupling model in 
the early universe, but deviates from all the other models at recent epochs due to the DE bounce. 



2.2 Linear perturbations 

For all the cosmological models considered in the CoDECS suite 
we need to compute, besides the specific evolution of relevant back- 
ground quantities like the Hubble function and the CDM mass vari- 
ation, as described above, also the evolution of linear density per- 
turbations in the Newtonian limit of sub-horizon scales. The growth 
factor D+(z) of each model will be in fact necessary in order to 
properly set up the initial conditions for the N-body runs described 
in Section [3] Although some fitting functions for the growth factor 
in the context of cDE models have been proposed in the literature 
Pi Porto & Amendola|2008l|Baldi et al.|2010||di Porto et al.|2012) 
these are valid only for the standard case of a constant coupling and 
for monotonic self-interaction potentials, and even in these sim- 
plified cases have a limited range of validity in terms of baryonic 
fraction and coupling strength. For the precision level required by 
the present work it is therefore still necessary to rely on direct nu- 
merical computations of the growth factor for each specific model 
under investigation. Differently from most of the previous nume- 
rical works involving cDE scenarios (see e.g. Maccio et al. 12004) 
|Baldi et aL|2"oTol|Baldi|20TTbl|Li & Barrow|2011| l, all the CoDECS 
runs assume a common normalization of the amplitude of density 



perturbations at zcmb, consistent with the latest constraints on the 
scalar perturbations amplitude A s from the WMAP7 results ( |Ko-| 
matsu et al. 201 1 ). All the models included in the CoDECS suite will 
therefore have a normalization of linear perturbations fully compa- 
tible with observational constraints coming from CMB data alone. 

In order to compute the growth factor D+ (z) for all our mod- 
els, we need to solve the linear perturbations equations for baryon 
and CDM density fluctiations in the presence of a DE coupling. 
These have been derived in the literature by several authors (see 
e.g. |Amendola|2004"l|Pettorino & Baccigalupi|2008||Baldi|2011b| >, 
and in Fourier space and cosmic time t take the form: 



8 C = -2H 



5 b = -2H5 b + AnG [p b S b + Pc 8. 



S c + 4ttG [p b S b + p c S c T c ] , (10) 



(11) 



where for simplicity we have omitted the <j> dependence of the cou- 
pling /3 C (4>)- In Eqs. jlO|ll[ l the overdensity of CDM and baryons 
is defined as 5 c ,b = Sp c ,b,/pc,b- The factor T c is defined as 

r c = i + p 2 c (t) (12) 

and includes the effect of the fifth-force mediated by the DE scalar 
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field for CDM density perturbations. The second term in the first 
square bracket at the right-hand-side of Eq. l |10| >, instead, rep- 
resents the so-called "extra friction" arising for CDM perturba- 
tions as a consequence of momentum conservation. The validity of 
Eqs. ( | 1 0| 1 1 1 > for the scales of interest in the CoDECS suite has been 
already verified in the literature both for the exponential ( |Amen-| 
|dola|2004||Baldi|2011bF and for the SUGRA bouncing cDE models 
l |Baldi|2012| >~ 

For each of the models under investigation we therefore solve 
numerically Eqs. < | 1 0| 1 1} along the corresponding background so- 
lution assuming as a boundary condition for the integration that the 
ratio of baryonic to CDM perturbations at zcmb takes the value 
Sb/5c ~ 3.0 x 10 -3 , as computed (for k > 1 h / Mpc) by runnin 



the publicly available Boltzmann code CAMErllLewis et al. 



2000 



for a ACDM cosmology with the WMAP7 parameters adopted in 
the present study (and listed in Table [TJ. We then define the to- 
tal growth factor D + (z) as a combination of the baryon and CDM 
perturbations amplitude weighted by the relative density fraction of 
the two components, as: 



D+(z) 



n c (z) + n b ( z ) 



[n c {z)s c (z) + n b (z)s b ( z )] . (13) 



The results of our integration are illustrated in Fig. [2] where 
we show the evolution of the normalized growth factor D+ / a and 
of its ratio to the ACDM case in the left and right panels, respec- 
tively. As one can see from the figures, all the models start by con- 
struction with the same amplitude of density perturbations at CMB, 
shown as a purple squared point in the plots, and then evolve with 
different growth rates until z = 0, thereby reaching significantly 
different normalizations of the linear perturbations amplitude at the 
present time. This implies a different value of erg at z = among 
the different cosmologies, as specified in Table [2] and as clearly 
visible in the figures, while all models are consistent with the same 
WMAP7 constraints at zqmb • 

Particularly interesting are the cases of the exponential cou- 
pling model EXP008e3 (orange lines) and of the bouncing cDE 
model (red lines). The former, as clearly shown in the right panel of 
Fig. [2] features a mild enhancement of linear perturbations growth 
for most of the expansion history of the universe, followed by a 
much more significant increase of the growth rate with respect to 
ACDM at low redshifts, when the coupling exponentially grows to- 
wards its present large value f3o = 0.4. The latter model, instead, 
as discussed in full detail by Baldi (2012), follows a kind of op- 
posite trend: after an extended initial phase of strongly enhanced 
growth, the evolution of density perturbations slows down dramati- 
cally as a consequence of the inversion of the scalar field motion at 
Zinv = 6.8 such that the standard ACDM amplitude is recovered at 
z — 0. Therefore, the bouncing cDE model SUGRA003 considered 
in this work has the same normalization of linear density perturba- 
tions as ACDM both at zcmb and at z = 0, while featuring a peak 
in its deviation from ACDM around z ~ 10. This peculiar evolu- 
tion is clearly visible by eye in both panels of Fig.[2]just by follow- 
ing the different paths of the black (ACDM) and red (SUGRA003) 
lines, and is expected to have significant consequences on structure 
formation at different cosmic epochs (Bald i|20l"2| >. 



3 THE CODECS SIMULATIONS SUITE 

As already mentioned above, the main purpose of the present pa- 
per is to present and discuss a vast numerical program for inter- 
acting DE cosmologies that has been recently carried out by the 
author and that is now made available to the scientific community 
for public use. We refer to this numerical enterprise as the COupled 
Dark Energy Cosmological Simulations (CoDECS) project, and in 
the present Section we will discuss the main features of the nume- 
rical implementation adopted for the different sets of runs carried 
out in the context of this project. 



2 www.camb.info 



3.1 Initial Conditions 

Initial conditions for N-body simulations are generated by dis- 
placing particles from a homogeneous distribution in order to set 
up a random-phase realization of the linear matter power spec- 
trum of the cosmological model under consideration according 
to Zel'dovich approximation ( Zel'do vichl 1970) . The particles dis- 
placements are then rescaled to the desired amplitude of density 
perturbation at some high redshift Zi when all perturbation modes 
included in the simulation box are still evolving linearly. This red- 
shift is then taken as the starting redshift of the simulation, and the 
corresponding particle distribution as the initial conditions for the 
N-body run. 

For all the simulations considered in the CoDECS project we 
generate initial conditions at Zj = 99 by displacing particles from 
a homogeneous glass distribution (White 1994 , Baugh et al. 1995 1. 
The initial power spectrum corresponds to the one obtained from 
CAMB at z — 99 for a ACDM cosmology with the cosmolog- 
ical parameters of Table [T| and includes baryonic acoustic oscil- 
lations (BAO). The same power spectrum is used for both CDM 
and baryon particles, which implies discarding possible effects due 
to a residual offset between the amplitude of perturbations in the 
two fluids at z — 99 (see e.g. Naoz et al. 2010), and for all the 
different cosmologies under consideration, thereby not consider- 
ing possible distortions of the transfer function due to early effects 
of the coupling. This choice is made in order to allow for a di- 
rect comparison of the impact that cDE cosmologies have on the 
growth of structures between zcmb and the present. Furthermore, 
the discarded effects are expected to have a much weaker impact 
on structure formation processes as compared to the DE interac- 
tions under investigation (Baldi e t al.|2010) . The particles displace- 
ments obtained for each simulation as a random-phase realization 
of the same initial power spectrum are then rescaled to the same 
amplitude at zqmb, corresponding to the WMAP7 normalization 
of scalar perturbations A s = 2.42 x 10 -9 . As we also employ the 
very same random seed for all the simulations, the resulting particle 
displacements are exactly the same for all the cosmological mod- 
els at zcmb and the initial conditions for each individual run are 
then generated by rescaling forward the displacements amplitudes 
to Zi = 99 with the correct growth factor D+ (z) (see Fig. \2.2\ 
computed for each cosmological model as described in Sec. [2] 

Once particles positions have been assigned, velocities are 
computed in Fourier space according to the linear perturbation the- 
ory relation v(k, a) oc f(a)8(k, a), where the growth rate function 
f(a) = din D+/ din a can be derived for each model from the 
corresponding computed growth factor. 
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Figure 2. Left: The growth factor (2) divided by the scale factor a for all the models included in the CoDECS project, normalized at the redshift of the 
last scattering surface, zqmb ~ 1100. Right: The ratio of the growth factor of each cDE model to the standard ACDM case. All models start with the same 
amplitude as ACDM at zqmb but significantly deviate from the standard evolution at later times. All the models show a monotonic increase of the deviation 
from ACDM, with the exception of the bouncing cDE model SUGRA003 that after an extended period of strongly enhanced growth inverts the trend in 
correspondence to the DE bounce, and recovers the ACDM amplitude at 2 = (see Baldi 201 lb 1. 



3.2 The N-body Simulations 

The CoDECS suite includes at the present time the six different 
cosmological models listed in Table [2] that have been extensively 
described in Sec. [2] For all these models we have run so far 
two different sets of cosmological N-body simulations, named 
L-CoDECS and H-CoDECS. Both sets of simulations consist of a 
cosmological volume with periodic boundary conditions filled with 
an equal number of CDM and baryonic particles, but differ from 
each other for the scale and for the physical processes included in 
the runs. All simulations have been carried out with the modified 
version by Baldi et al. 1201 0} of the wi dely used parallel Tree-PM 
N-body code GADGET l |Sprmgel|2005| >, specifically developed to 
include all the additional physical effects that characterize cDE 
models (see |Baldi et al.||2010] for a detailed description of the 
implementation of cDE in the code). 

The L-CoDECS simulations have a box size of 1 comoving 
Gpc/h aside and include 1024 3 CDM and baryon particles for a 
total particle number of 2 x 1024 3 w 2 x 10 9 . The mass resolution 
at z = for this set of simulations is m c = 5.84 x 10 10 M //i 
for CDM and m b = 1.17 x 10 10 M Q /h for baryons. Despite the 
presence of baryonic particles these simulations do not include hy- 
drodynamics and are therefore purely collisionless N-body runs. 
The inclusion of baryonic particles is necessary in order to real- 
istically follow the growth of structures in the context of specific 
cosmological scenarios (as the cDE models under discussion here) 
where baryons and CDM do not obey the same dynamical equa- 
tions. In other words, if the fraction of uncoupled matter given by 
the baryons was omitted, the effect of the coupling on structure 
formation would be significantly overestimated in the simulations 
leading to largely incorrect results. On the other hand, the different 
dynamics followed by CDM and baryonic particles imprints chara- 
cteristic features on the relative distribution of the two fluids which 
might have interesting observational implications, and therefore the 
inclusion of collisionless baryonic particles in the runs offers also 
the possibility to estimate such effects, as we will discuss in detail 
in Sec.l4~2l 

In Fig. [5] we show the CDM density distribution at z = in 
a slice of 1000 x 250 Mpc/ft for the standard ACDM cosmology 



(middle) and for two of our cDE models, namely the most extreme 
model with exponential potential and constant coupling EXP003 
(top) and the bouncing cDE model SUGRA003 (bottom). The 
slice is 30 Mpc/h thick and is centered on the most massive halo 
identified for the ACDM run. Although the main pattern of the 
large-scale cosmic web looks very similar in the three models 
(as a consequence of the identical random phases adopted for 
the generation of the initial conditions in each simulation), some 
differences can be identified even by eye when comparing the 
density contrast of the most concentrated objects as well as of the 
largest voids. The EXP003 model appears clearly more structured 
and with a more pronounced density contrast between halos (in 
yellow) and void regions (in black) as compared to ACDM. Fur- 
thermore, by carefully looking at the central part of the two slices, 
it is also possible to notice how some of the satellite halos that are 
falling along the filaments onto the main central structure in the 
ACDM cosmology have already merged into bigger substructures 
in the EXP003 case. The SUGRA003 model shows much weaker 
signatures with respect to ACDM at z = 0. This is expected 
since the model has the same as normalization as ACDM at the 
present time. A more quantitative estimate of the differences in 
the abundance and distribution of massive halos is provided by 
the Halo Mass Function extracted from the L-CoDECS runs that is 
briefly sketched in Sec. |4.3| and in Figure[7j and discussed in much 
finer details in |Baldi|p012( and |Cui et al.|j2012) . 

The H-CoDECS simulations are instead adiabatic hydrodynam- 
ical simulations in a cosmological box of 80 Mpc/h aside filled 
with 512 3 CDM and gas particles, where the gas particles repre- 
sent the uncoupled baryonic fraction. The mass resolution at 2 = 
for this set of simulations is m c = 2.39 x 10 s Mq/H for CDM 
and rut — 4.78 x 10 7 M jh for baryons. Adiabatic hydrodynam- 
ical forces on the gas particles are computed by means of the en- 
tropy conserving formulation of Smoothed Particle Hydrodynamics 
(SPH, Springel & Hernquist 2002) and other radiative processes 
such as gas cooling, star formation, or feedback mechanisms are 
not included in the simulations. 

In Fig. [4] we show the formation of a massive galaxy clus- 
ter in the same three cosmological models displayed in Fig. [5] by 
plotting the gas density contrast in a slice of 40 x 40 Mpc/h cen- 
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Figure 3. The CDM density distribution in a slice with size 1000 X 250 Mpc//i and thickness 30 Mpc//i as extracted from the L-CoDECS simulations of 
a few selected models. The middle slice shows the case of the standard ACDM cosmology, while the top slice is taken from the EXP003 simulation and the 
bottom slice from the bouncing cDE model SUGRA003. While the latter model shows basically no difference with respect to ACDM at z = 0, due to the 
very similar value of erg for the two models, clear differences in the overall density contrast and in the distribution of individual halos can be identified by eye 
for the EXP003 cosmology. (A higher resolution version of this figure is available online through the CoDECS website, see the Appendix) 



tered in the position of the most massive halo identified at z — 
in the ACDM simulation of the H-CoDECS suite, for three differ- 
ent redshifts z = 0, 1, and 2. In this case, since the figure shows 
the redshift evolution of the three models, it is possible to identify 
the differences that characterize the growth of structures in these 
selected scenarios. Also this comparison shows the same features 
already described above for the comparison between ACDM and 
EXP003. By comparing the two cosmologies at z = 0, in fact, it 
clearly appears how the latter has a more evolved structure of the 
cosmic web around the central massive halo and shows a larger 
density contrast between the highly nonlinear overdense regions 
and the voids. It also appears much more clearly than before, from 
these plots, the lack of small structures in the EXP003 universe 
with respect to ACDM, due to the faster merging processes occur- 
ring in the cDE cosmology. Despite this, the relative abundance of 
substructures as a function of their fractional mass with respect to 
their host main halo is not significantly affected by the coupling, 
as will appear from the Subhalo Mass Function displayed in Fig. [5] 
below. 



The higher amplitude of density perturbations in the EXP003 
model is visible also at higher redshifts, especially at z = 1, where 
the difference in the structure of the cosmic web around the central 
halo between EXP003 and ACDM is particularly evident, as well as 
the different merging history of the two models. Particularly inter- 
esting is in this case also the comparison between ACDM and the 
bouncing cDE model SUGRA003. At z = 0, as already discussed 
above, the two models show much weaker differences with respect 
to the case of the EXP003 cosmology, due to the comparable am- 
plitude of their linear density perturbations. However, at high red- 
shifts (e.g. at z — 2) the SUGRA003 model appears clearly more 
evolved as compared to ACDM: the main filamentary structure and 
the high density peaks appear more luminous than the correspond- 
ing ACDM case, and the main halos are more extended due to the 
merging of the surrounding satellite structures. This shows the main 
characteristic feature of the bouncing cDE scenario: a larger ampli- 
tude of density perturbations (and a correspondingly higher number 
of massive clusters) with respect to ACDM at high redshifts, fol- 
lowed by a reconnection of the two models that show very little 
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Figure 4. The gas density distribution during the formation process of a massive galaxy cluster as extracted from the H-CoDECS runs for the same three 
models shown in Fig. [3] Also in this case, differences in the overall density contrast and in the distribution of individual lumps are visible by eye when 
comparing the standard ACDM cosmology and the EXP003 cDE model at z = 0. However, in this case the redshift evolution shown in the figure allows 
to identify differences also between ACDM and the bouncing cDE model SUGRA003 at higher redshifts, where the latter model appears more evolved and 
shows a more pronounced density contrast as compared to the standard cosmology. (A higher resolution version of this figure is available online through the 
CoDECS website, see the Appendix) 

differences in their large-scale structures at the present time ( |Baldi| 4 BASIC RESULTS 

|2012l >. 

We describe here some basic results arising from the analysis of 
the CoDECS simulations, focusing on the nonlinear matter power 
High-resolution versions of the images displayed in spectrum, on the effects of the coupling on the relative evolution of 
Figs. [5] and [4] are available online through the CoDECS web- CDM and baryon density perturbations, and on the halo and sub- 
site (see the Appendix). halo mass functions. As the simulations outputs are now publicly 
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available, more detailed and extended analysis might be carried out 
in the near future also by other independent groups. This is the main 
aim of the CoDECS public release discussed in the Appendix to this 
paper. 

4.1 The nonlinear matter power spectrum 

For the L-CoDECS runs we compute the density power spectrum 
of the different components P b (k), P c {k), Pm(&) (for the CDM, 
baryon, and the total matter components, respectively) by deter- 
mining the individual density fields on a grid with the same size 
of the PM grid (i.e. 1024^ nodes) through a Cloud-in-Cell mass as- 
signment. This procedure provides a determination of the nonlinear 
matter power spectrum up to the Nyquist frequency of the PM grid, 
which for the L-CoDECS runs corresponds to &N y = nN/L as 
3.2/i/Mpc. Beyond this limiting wave number we estimate the 
power spectrum by means of the folding method of Jenkins et al. 
( 1998 ); Col ombi et al.|p009] l, and we smoothly interpolate the two 
estimations around &Ny in order to obtain a single determination of 
the matter power spectrum at different redshifts for all the simula- 
tions of the sample. The small scale power spectrum is then trun- 
cated at the frequency where the shot noise reaches 10% of the 
measured power. 

With the nonlinear power spectra determined as discussed 
above, we can investigate the effects that the coupling induces on 
the growth of structures and highlight how differently these effects 
characterize the linear and the nonlinear regimes of structure for- 
mation. We should stress once more that the CoDECS runs differ 
from all the previous N-body simulations of interacting DE models 
as they feature a common normalization of density perturbations at 
CMB, while previous studi es jMaccio et al.|2004"l|Baldi et al.|2010[ 



|Li & B arrow 2011; Baldi 2011b) have almost always assumed a 
common normalization of linear perturbations at z = 0. Our re- 
sults therefore represent the first realistic estimation of the effects 
of cDE models for a growth history of density perturbations fully 
compatible with CMB data. 

As a consequence of the common normalization of perturba- 
tions at CMB and of the different growth factors determined by the 
modified expansion history and by the fifth-force and extra friction 
terms in Eq.[l0] the various cDE models will have in general a dif- 
ferent amplitude of density perturbations at any redshift z < zcmb 
and in particular at the present time z — 0. This will imply a dif- 
ferent value of erg from model to model, as summarized in Table|2] 
with the only exception of the SUGRA003 bouncing cDE model 
that has been constructed in order to have the same linear normal- 
ization as ACDM both at zcmb and at the present epoch. 

Furthermore, the coupling between DE and CDM has very 
different effects on the growth of structures in the linear and in 
the nonlinear regimes. In particular, as it has been shown in pre- 
vious studies (see e.g. |Baldi|[201 la| for a detailed discussion of 
this effect), the extra friction term of Eq.[l0]for the case of stan- 
dard cDE models (as the exponential models with constant cou- 
pling discussed in this work) accelerates CDM particles along their 
direction of motion, and therefore acts in the direction of enhancing 
the growth of structures in the linear regime, while in the nonlin- 
ear regime, due to the onset of virialization processes and the ap- 
pearance of non-radial velocities in the local peculiar flow of parti- 
cles, it slows down the contraction of collapsed structures thereby 
slowing down structure formation processes. This determines a re- 
duction of the small scale nonlinear power in cDE models with 
respect to the expected power based on their large-scale linear nor- 
malization. Such effect is completely inverted for the case of the 



SUGRA bouncing cDE model, where the sign of the extra-friction 
term changes in correspondence to the "bounce" of the DE scalar 
field at Zi nv ~ 6.8. For this model, in fact, since c/>/3 c (cA) < for 
z < z- mv ~ 6.8, the friction term acts as a proper friction during 
the late stages of structure formation, thereby slowing down parti- 
cles along the direction of their motion. This implies that the fric- 
tion term alleviates the effect of the fifth-force by slightly slowing 
down the growth of structures in the linear regime at low redshifts, 
while it favors the contraction of collapsed objects in the nonlinear 
regime. 

This opposite behavior of the extra-friction term in the linear 
and nonlinear regimes determines specific types of distortion of the 
matter power spectrum that allow to break the degeneracy of the 
effect of the coupling at the linear level with other cosmological 
parameters - as e.g. erg - as shown in Fig. [5] In each of the plots of 
Fig. [5] in fact, we display with solid lines the ratio of the nonlinear 
matter power spectrum of each cDE model to the fiducial ACDM 
cosmology with erg = 0.809 as extracted from the L-CoDECS 
simulations. As a useful comparison, the dotted line shows the 
ratio of the nonlinear matter power spectra computed with the 
Boltzmann code CAMB by means of the HALOFIT fitting routine 
(Smith et al. 2003) for a ACDM model with a erg normalization 
corresponding to the one of the respective cDE model and for the 
fiducial ACDM cosmology. In other words, the solid lines show 
the ratio P sim (cDE)/P sim (ACDM) while the dotted lines show 
the ratio Pcamb (cDE as)/ Pcamb (fiducial <t s ). 

The three upper plots of Fig.|5]display these ratios for the three 
standard cDE models EXP001, EXP002, and EXP003 at different 
redshifts. For very small couplings (as e.g. for EXP001, cyan line) 
the effect on the matter power spectrum at z = appears to be 
fully degenerate with as, as the two ratios look indistinguishable 
at all scales. However, for larger couplings (EXP002, green line) a 
suppression of power in the cDE model with respect to a ACDM re- 
alization with identical erg is clearly visible at small scales. The ef- 
fect becomes more prominent for the strongest coupling (EXP003, 
blue lines), where a significant deviation between the cDE model 
and a ACDM cosmology with as = 0.967 appears at z — for 
k > 0.3 /i/Mpc. Therefore, even at z — the nonlinear scales of 
the matter power spectrum allow to distinguish between a cDE sce- 
nario and a ACDM model with the same large-scale amplitude of 
linear perturbations, removing the apparent degeneracy that arises 
when only linear scales are considered. This difference might be 
beyond the discriminating power of present observational probes 
for very small couplings, but represents in any case a distinctive 
signature of cDE scenarios. 

Furthermore, by looking at the redshift evolution of the 
large-scale matter power in the three upper plots of Fig. [5] it is 
clearly possible to distinguish, at least in principle, a cDE model 
from a ACDM cosmology with high erg by the different evolution 
of the amplitude of density perturbations associated to the different 
linear growth factors displayed in Fig. [2] All the standard cDE 
models EXP001-EXP003, in fact, show a lower amplitude of 
the linear matter power at high redshift as compared to what ex- 
pected for a ACDM cosmology with the same value of erg at z = 0. 

The two lower plots of Fig.[5]show the same quantities for the 
exponential coupling model EXP008e3 and for the bouncing cDE 
model SUGRA003. Interestingly, while the former shows a simi- 
lar behavior as the standard cDE models with constant coupling, 
with a large-scale amplitude at z = fully consistent with the 
linear growth expected for this model, and with a small-scale sup- 
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Figure 5. The ratio of the nonlinear matter power spectrum to the standard ACDM cosmology for all the cDE models under discussion, as extracted from 
the L-CoDECS runs (solid), and the ratio of the nonlinear matter power spectrum for a ACDM model with the same erg as the respective cDE model and the 
fiducial ACDM cosmology with erg = 0.809, as computed by means of the HALOFIT routine of the public code CAMB (dotted). Although at linear scales 
the two quantities are completely degenerate at 2 = 0, the different redshift evolution of the linear power spectrum amplitude allows to break the degeneracy. 
Furthermore, even at 2 = 0, the nonlinear regime shows a clear difference between the effects of cDE and of a high value of erg due to the distortion determined 
by the DE coupling on the matter power spectrum, be it either a suppression (as for constant and variable coupling models) or an enhancement (as for bouncing 
cDE). 



pression of power due to the nonlinear effect of the friction term, 
the latter model is found to have the opposite trend both concern- 
ing the redshift dependence of the large-scale perturbations ampli- 
tude and the small-scale distortion of the nonlinear matter power. 
In particular, the linear amplitude of density perturbations is fully 
consistent with the fiducial ACDM cosmology at the present time, 
while a significantly higher normalization with respect to the fidu- 
cial model clearly appears at high redshifts. This offers a direct way 
to constrain this scenario and to account for a possible anomalous 
clustering at high z (Baldi 2012) with the next generation of large 
galaxy surveys. At the nonlinear level, instead, the model shows a 
strong increase of the small scale power at low redshifts in a range 
of scales (k ~ 1 — 10/i/Mpc) potentially accessible to accurate 
weak lensing measurements. 

Besides its peculiar impact on the redshift evolution of the 
cluster mass function (see Baldi 2012, for a study of the expected 
cluster number counts), that remains the main motivation behind 
the whole bouncing cDE scenario, these specific signatures identi- 
fied for the first time by the present study clearly provide a handle 
to observationally test this class of cosmologies. 



4.2 The gravitational bias 

A second characteristic feature of cDE models widely discussed in 
the literature (see e.g. Ma ccio et al.|2 004, Maini m~& Bonometto| 



|2006|[Baldi et al.|2010) arises as a consequence of the different dy- 
namic equations that CDM and baryonic particles obey due to their 
different couplings to the DE scalar field. While baryons are al- 
ways uncoupled and therefore follow standard Newtonian dynam- 
ics, CDM particles are affected by additional forces (the fifth-force 
and the extra friction) due to the interaction with the DE scalar 
field, as shown in Eqs. < | 1 0| 1 1[ >. This different dynamical evolution 
of the two matter components of the universe determines an off- 
set in the amplitude of density perturbations of the two fluids, that 
starting from the same amplitude at high redshifts progressively de- 
velop a "gravitational bias", which has been shown to have a very 
significant impact on the expected baryon fraction of bound objects 
as galactic or cluster sized halos (Bal di et al.|20 10 Baldi 2 01 lb| >. 
More specifically, a lower baryon fraction is genetically expected 
in cDE models as compared to ACDM as a consequence of the seg- 
regation of baryons (first discussed by Mainini & Bonometto |2006| > 
in the outskirts of collapsed structures determined by the gravita- 
tional bias. The redshift and scale dependence of this effect can be 
tested by computing the ratio of the nonlinear matter power spectra 
of the baryonic and CDM components, Pb(k)/ P c (k). In Fig.[6]we 
plot this quantity, normalized to the ACDM case, for all the cDE 
models considered in the present work and at three different red- 
shifts. As one can see from the figure, the gravitational bias shows 
a very characteristic scale dependence: while in the linear regime 
represented by the largest scales a weak and scale-independent off- 
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Figure 6. The ratio of the baryon to CDM power spectra as a function of 
scale and at different redshift: 2 = 1 (top), z = 0.56 (middle), z = 
(bottom). The nonlinear enhancement of the bias clearly develops as time 
passes and strongly reduces the baryon fraction of overdense regions at 
small scales, while at the largest scales the effect never exceeds a few per- 
cent. The growing coupling model EXP008e3 and the bouncing cDE model 
SUGRA003 show the largest nonlinear bias, which nevertheless is strongly 
enhanced for fundamentally different reasons in the two cases, as explained 
in the text. 



set appears for all the models, with a maximum deviation from the 
standard ACDM value of a few percent, the effect becomes much 
larger and strongly scale-dependent when moving to progressively 
more nonlinear scales, reaching at z = a suppression of the bary- 
onic power with respect to the CDM power (at k ~ 10 /i/Mpc) of 
about 1%, 5%, and 11% for the standard cDE models with constant 
coupling EXPOOI (cyan), EXP002 (green), and EXP003 (blue), re- 
spectively. 

Particularly interesting are again, however, the cases of the ex- 
ponential coupling model EXP008e3 (orange) and of the bouncing 
cDE model SUGRA003 (red). These show a much steeper scale 
dependence in the mildly-nonlinear and nonlinear ranges of scales 
as compared to the standard cDE models with constant coupling, 
reaching a suppression of the baryonic power with respect to the 
CDM power of ~ 30% and ~ 15%, respectively, while at large 
scales the effect has a comparable amplitude to the standard cDE 
models EXPOOI -EXP003. Interestingly, this enhancement of the 
gravitational bias at small scales that characterizes the EXP008e3 
and the SUGRA003 models has a different origin in the two cases. 

The exponential coupling model EXP008e3, in fact, features 
a sudden enhancement of the growth of CDM perturbations due 
to the steep increase of the coupling /3 c (<f>) at low redshifts. This 
determines a fast increase of the effective gravitational constant 
G cff = GT c (<f)) (see Eqs. that applies to CDM particles 
only, thereby inducing a fast collapse of CDM overdensities with 
respect to the corresponding baryonic perturbations. The latter are 
in fact governed by the standard gravitational constant G and there- 
fore respond with a longer dynamical timescale to the rapidly grow- 
ing gravitational potential induced by the CDM collapse. As a con- 
sequence, the growth of baryon fluctuations is always delayed with 
respect to the corresponding CDM perturbations, giving rise to 
the large bias observed at small scales. This interpretation is con- 
firmed by following the redshift evolution of the effect: the scale 
dependence of the bias and the amplitude of the small-scale offset 
strongly increase between 2 = 1 and z — 0, in correspondence 
to the steep increase of the coupling function for this model. This 
effect was already investigated by Baldi (201 lb I for a wide range 
of growing coupling scenarios with a common normalization of 
linear perturbations at z = 0. We therefore confirm here that the 
same effect appears also for models with a common normalization 
of linear perturbations at CMB, although with a significantly larger 
amplitude than for the former case. 

The bouncing cDE model SUGRA003, instead, does not show 
a very strong redshift evolution of the scale dependence of the grav- 
itational bias. In this case, in fact, the origin of the effect is substan- 
tially different than for the exponential coupling model. As dis- 
cussed in full detail by Baldi (2012), the inversion of the direction 
of motion of the DE scalar field in the bouncing cDE scenario de- 
termines a change of sign of the quantity 4>/3 c , that for the specific 
model considered here goes from positive to negative in correspon- 
dence to the DE "bounce" at Z[ nv ~ 6.8. This implies a change of 
sign of the extra friction term of Eq. < | 1 0[ > that behaves as a drag 
term (i.e. accelerating particles along their direction of motion) be- 
fore the bounce and as a proper friction (i.e. decelerating particles) 
after it. In the low-redshift universe, then, the SUGRA003 model is 
characterized by an extra friction term that decelerates CDM parti- 
cles thereby favoring the collapse of CDM perturbations by remov- 
ing angular momentum only from the CDM component. This de- 
termines a contraction of CDM perturbations that grow faster than 
their baryonic counterparts, giving rise to the observed enhance- 
ment of the gravitational bias with respect to standard cDE scenar- 
ios. 
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The level of baryon depletion of density perturbations in the 
different cDE models will be reflected in the baryonic fraction of 
massive halos (see Baldi et al.|201Q[ [Baldi 20 lib) and could pro- 
vide a further direct observational signature of cDE scenarios. It 
is important to notice that the plots are taken from the L-CoDECS 
suite which does not include any hydrodynamical force. The off- 
set between the two components is therefore of pure gravitational 
origin, due to the effective modification of the gravitational dynam- 
ics of CDM particles in cDE models. Additional physical processes 
as hydrodynamical pressure, gas cooling, star formation, and feed- 
back mechanisms are therefore expected to modify this picture in a 
substantial way (see e.g. Duffy et al. 2010; van Daaler Tet al.|201 1| 
|Semboloni et al.||2011| >. However, it should be stressed here that 
the gravitational bias described in the present section is a purely 
dynamical mechanism that starts to develop well before any of 
such effects can take place, thereby reducing the baryonic reservoir 
available for astrophysical processes. Detailed high-resolution sim- 
ulations of cDE models including all these additional ingredients 
will therefore be necessary in order to provide a reliable quantita- 
tive determination of the impact of the gravitational bias on observ- 
able properties of bound halos. Nevertheless, the simple estimate 
discussed in the present Section already provides a clear descrip- 
tion of the redshift and scale dependence of the expected bias. 

4.3 The Halo Mass Function 

The enhanced growth of density perturbations and the modified 
nonlinear dynamics of CDM particles have a significant impact on 
the expected number of halos as a function of redshift and mass, en- 
coded by the Halo Mass Function (HMF) N(> M, z). The impact 
of cDE cosmologies on the HMF has been studied in the past by 
several authors (s ee e.g. |Macci6 et al.|2004||Baldi et al.|2010||Baldi| 
& Pettorino 201 1 , Li & Barrow 201 1 and references therein), while 
a specific investigation of the expected number counts of massive 
clusters at high redshift in the context of the bouncing cDE sce- 
nario SUGRA003 discussed in this work has already been carried 
out by Baldi (2012). Finally, a detailed and comprehensive analy- 
sis of the HMF in cDE models based on the CoDECS simulations 
has been recently presented by |Cui et al.| ( |2012} , resorting both on 
the Friends-of-Friends (FoF) and on a series of Spherical Overden- 
sity halo samples for all the cosmological models considered in the 
present work. 

Therefore, we refer the interested reader to the latter work for 
a thorough discussion on the features of the HMF in cDE scenar- 
ios, and we present here only a simple comparison of the expected 
number counts of massive halos in the different cosmologies con- 
sidered in the CoDECS project based on the FoF sample. This has 
been constructed by running a FoF algorithm with linking length 
I = 0.2 X d, where d is the mean interparticle separation (see 
Appendix [A} on the CDM particles as primary tracers of the lo- 
cal mass density, and then attaching baryonic particles to the FoF 
group of their nearest neighbor. In Figure [7] we show the HMF ex- 
tracted from the L-CoDECS runs which has been computed using 
the FoF catalogs at high redshift (z — 2.5, upper panel) and at the 
present epoch (z = 0, lower panel) for all the cosmological mod- 
els included in the CoDECS project. The bottom plot in each panel 
displays the ratio of the number counts with respect to the fiducial 
ACDM cosmology. 

As one can see from the figures, all the cDE models pre- 
dict a significant enhancement in the number of halos at high red- 
shift. The effect is visible over the whole mass range covered by 
our FoF sample, but shows a clear mass dependence with a much 
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Figure 7. The Halo Mass Function of the different cDE models included 
in the CoDECS project at z = 2.5 (upper panel) and at z = (lower 
panel) as extracted from the halo catalogs of the L-CoDECS runs based 
on a Friend-of-Friend identification algorithm. The enhancement of the ex- 
pected number counts of massive halos is evident in both plots, with a clear 
mass dependence of the effect. All the models show a significant increase 
of the number of very massive objects at the present epoch except for the 
case of the "Bouncing" cDE model SUGRA003 (see also the more detailed 
analysisof the Halo Mass Function in cDE models presented in Baldi 2012 
|Cui et al.|20T2> 



stronger enhancement at large masses, reaching a factor of ~ 20 at 
M w 1.0 x 10 14 M Q /ft for the most extreme scenario EXP003. 
The other models have a weaker - although still significant - impact 
on the expected number counts of massive halos, and it is particu- 
larly interesting to notice the enhancement of a factor of ~ 4 at 
M » 1.0 x 10 14 M©/ft of the bouncing cDE model SUGRA003 
(see also |Baldi|2012) . 

At low redshift the effect becomes more strongly mass depen- 
dent and moves towards larger masses. The most extreme scenario 
shows an enhancement of a factor of ~ 15 at M w 3.0 x 10 15 
Mg/A, which reduces to a factor of ~ 2 — 5 for the other mod- 
els. The only exception is provided by the bouncing cDE model 
SUGRA003, that shows basically no difference with respect to 
ACDM at z — (see again the detailed discussion in Baldi [2012) . 

These results show some of the possible applications of the 
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CoDECS public halo catalogs that are introduced in the Appendix. 
The enhancement of the expected number of massive clusters as a 
function of redshift is one of the distinctive features that character- 
ize cDE scenarios and might provide a direct way to constrain the 
parameters of the model. 

4.4 The subhalo mass function 

As a further example of the possible applications of the CoDECS 
public catalogs we present here the subhalo mass function, defined 
as the number of substructures with a given fractional mass with 
respect to the virial mass of their host main halo (M su b/M2oo) 
as a function of the fractional mass itself. We compute this quan- 
tity using the subhalo catalogs obtained for all the simulations of 
the CoDECS project by means of the SUBFIND algorithm ( Springel 
|et al.||200T) , and subsequently binning in logarithmic fractional 
mass bins the whole sample of substructures belonging to host ha- 
los with virial mass above a minimum mass threshold. In Figure [8] 
we show the subhalo mass function at z — as extracted from the 
subhalo catalogs of the H-CoDECS simulations. The bottom plot 
displays the ratio of the subhalo mass function with respect to the 
fiducial ACDM cosmology. The mass threshold for the main halos 
has been chosen as M m m = 5 x 10 12 Mg/ft. 

As one can see from the figure, the mass distribution of halo 
substructures is not significantly affected by the interaction be- 
tween DE and CDM, as no clear trend comparable to what observed 
for the HMF appears in the plot, besides some random scatter at 
large fractional masses due to the limited number of subhalos in the 
high-mass bins of the sample. This result shows for the first time 
how the coupling between DE and CDM does not affect the rela- 
tive abundance of halos of different masses but rather determines 
an overall shift in the mass of gravitationally bound structures dur- 
ing the whole assembly process of large halos. A more detailed 
analysis of the formation process of such halos and of the spatial 
distribution and dynamical state of subhalos, which goes beyond 
the scope of the present paper, would be required in order to detect 
possible footprints of DE interactions in the observable features of 
halo substructures. 



5 CONCLUSIONS 

We have described the set of cosmological models and the nume- 
rical setup adopted for a vast simulations program aimed at inves- 
tigating the nonlinear evolution of Coupled Dark Energy cosmolo- 
gies to an unprecedented level of statistical significance and nume- 
rical accuracy. We refer to such enterprise as the "CoDECS project", 
standing for COupled Dark Energy Cosmological Simulations. The 
CoDECS project includes different types of interacting Dark Energy 
scenarios and covers a wide range of scales by means of two sepa- 
rate suites of cosmological N-body simulations on large scales (L- 
CoDECS) and small scales (H-CoDECS) for each of the models un- 
der investigation, the latter including also an hydrodynamic treat- 
ment for the uncoupled fraction of baryonic particles. 

The range of models considered in the CoDECS suite includes 
constant coupling, variable coupling, and bouncing Coupled Dark 
Energy scenarios, and for each of these cosmologies we have de- 
scribed in full detail the background and linear perturbations evo- 
lution, derived by numerically integrating the system of coupled 
differential dynamic equations from the very early universe to 
the present time. Linear perturbations have been normalized to 
the same amplitude at the last scattering surface in accordance to 
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Figure 8. The Subhalo Mass Function at z = of all the cDE models in- 
cluded in the CoDECS project as extracted from the subhalo catalogs of the 
H-CoDECS runs based on the SUBFIND algorithm ^Springel et al.|2001fr . 
The plot shows the Subhalo Mass Function for host halos with virial mass 
M200 > 5.0 X 10 12 Mq//i. No clear trend appears in the relative abun- 
dance of subhalos as a function of the different realizations of DE interac- 
tions. 

the latest WMAP7 results. This allows to study, by means of the 
CoDECS simulations, the evolution of structure formation in the 
context of different Coupled Dark Energy scenarios that are all con- 
sistent with present CMB constraints. 

We have described in full detail the setup of initial conditions 
and the numerical features of our different N-body simulations, 
and provided a visual comparison of the structures forming in 
the context of a few selected models as extracted from both the 
L-CoDECS and H-CoDECS runs. Although the main geometrical 
structure of the cosmic web appears very similar in all the models, 
due to the identical phases adopted for the initial conditions 
generation, some differences in the overall density contrast and in 
the distribution and merging history of the smaller structures can 
be easily identified even by eye. In particular, for those models 
that due to the coupling end up with a significantly higher value of 
erg at z = 0, the difference in the density contrast between highly 
nonlinear lumps and void regions is particularly evident. 

We have then moved to a more quantitative analysis of the 
CoDECS simulations, by investigating the impact of the differ- 
ent cosmological models on the nonlinear matter power spectrum. 
More specifically, we have studied the degeneracy of the effects 
of the Dark Energy coupling with other cosmological parameters, 
in particular with erg, by comparing the nonlinear power spectra 
extracted from each of our simulations to the power spectra com- 
puted for a ACDM cosmology with the same erg of the considered 
Coupled Dark Energy model by means of the HALOFIT module 
implemented in the publicly available code CAMB. 

At linear scales, although a perfect degeneracy was found 
between Dark Energy coupling and <rg at z = 0, our results have 
shown for the first time how such degeneracy could be broken, 
at least in principle, by the redshift evolution of the linear power 
spectrum amplitude, that in Coupled Dark Energy decreases more 
rapidly with redshift as compared to a ACDM cosmology with 
the same erg at z = 0. Even more remarkably, we have shown 
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how the degeneracy can be easily broken, for sufficiently large 
couplings, even at z = by moving to nonlinear scales beyond 
k ~ 0.3 /i/Mpc: the nonlinear effect of the friction term associated 
to the Dark Energy interaction significantly suppresses the non- 
linear power for constant and variable coupling models, while it 
strongly increases it for bouncing Coupled Dark Energy scenarios, 
thereby providing a clear distinction between an interacting Dark 
Energy and a high-ag ACDM cosmology. 

We have also studied the effect of the Dark Energy coupling 
on the offset between the amplitude of density fluctuations in 
the CDM and baryonic components, the so called gravitational 
bias. Due to the different dynamics determined by the interaction 
with the Dark Energy scalar field, CDM perturbations are in 
fact expected to develop faster with respect to the corresponding 
baryonic fluctuations, that obey standard Newtonian dynamics. 
Our analysis has shown, confirming previous results, how the 
offset between the perturbations amplitude in the two compo- 
nents, that is limited to a few percent and does not show any 
scale dependence at linear scales, gets strongly enhanced in the 
nonlinear regime (k it 0.3 — 1.0 /i/Mpc), where it acquires also a 
strong scale dependence. Such effect has been shown in the past 
to have a direct impact on the predicted cluster baryon fraction, 
determining a lower baryonic content of massive CDM halos. Here 
we have highlighted the cases of growing coupling scenarios and 
of bouncing Coupled Dark Energy models, for which the effect 
of baryon segregation at small scales is particularly dramatic, 
thereby offering a further signature to observationally constrain 
such models. 

Finally, we have provided a few examples of the possible 
applications of the public CoDECS halo catalogs by computing 
the Halo Mass Function from the L-CoDECS runs at high redshift 
and at the present epoch and the Subhalo Mass Function from the 
H-CoDECS runs at z = 0. Our results have shown how the relative 
abundance of massive halos is enhanced in Coupled Dark Energy 
scenarios with respect to ACDM in a mass- and redshift-dependent 
way. Furthermore, we have shown for the first time that the mass 
distribution of substructures within large collapsed objects is not 
significantly different in Coupled Dark Energy models with respect 
to the standard ACDM cosmology. 

The main aim of the present paper is to provide a general 
description of the cosmological models and of the numerical setup 
of the CoDECS project, as a a reference for the public release of 
the simulations outputs. For this reason, we have also described in 
the Appendix of the paper the main features of the CoDECS public 
database, and the different types of data that are made directly 
available for public use. More detailed analysis of the CoDECS 
runs will be performed in future works. 

To conclude, we have performed the largest and most detailed 
numerical study of Coupled Dark Energy models to date, and pro- 
vided a general description of our numerical methods and findings. 
Our results broadly confirm previous outcomes but significantly 
extend the statistical significance of the results and the range of 
models covered by the simulations with respect to previous works, 
also including for the first time the recently proposed Bouncing 
Coupled Dark energy scenario. Furthermore, differently from all 
previous numerical studies of Coupled Dark Energy models, our 
CoDECS simulations are fully consistent with the latest CMB con- 



straints from WMAP7 concerning the normalization of both the 
background cosmological parameters and the amplitude of den- 
sity perturbations. The simulations outputs, including halo catalogs 
and nonlinear matter power spectra, are made publicly available 
through a web database and will serve for future analysis of inter- 
acting Dark Energy cosmologies. 
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APPENDIX A: THE CODECS PUBLIC DATABASE 

All the numerical outputs of the CoDECS project are available for 
public use. These include both the raw snapshots data of the indi- 
vidual simulations and some processed quantities as the nonlinear 
matter power spectra or the halo and sub-halo catalogs. The post- 
processed data are made directly accessible through a web interface 
at the URL: 

www.marcobaldi.it/research/CoDECS 

while the unformatted snapshots data, due to their large size, are 
stored on a separate server and must be requested directly to the au- 
thor. In the present section we will briefly describe which data are 
available and how post-processing has been performed, while we 
refer to the specific user's guide accessible at the CoDECS website 
mentioned above for a more specific description of how to interpret 
the file formats of the processed and raw data. The numerical fea- 
tures of the two different simulation suites included in the CoDECS 
project - the L-CoDECS and the H-CoDECS - are summarized in 
Table IaTI 

Al The CoDECS website 

The CoDECS website provides all the necessary information to ac- 
cess, download, and interpret the CoDECS data, including the terms 
of use and the relevant references. From the URL address given 
above, one can access a summary page where the main features of 
all the simulations of both the L-CoDECS and H-CoDECS suite are 
listed and from where a specific database for each simulation can 
be accessed. Each of the different simulations, in fact, has an indi- 
vidual web page containing all the directly available data for that 
specific run. All the individual pages have a standard format that is 
described in detail below. 

A2 Background and linear perturbations data 

For each specific model, a number of formatted files containing in- 
formation about the background and linear perturbations evolution 
are available. 

The background files: 

modeLCoDECS-background.dat 

provide the redshift evolution of the relative energy fraction of all 
the different components of the universe f2 c , f2&, Q r , f2ki n , f2 pot , 
where the latter two quantities represent the fractional energy 
density of the kinetic and potential contributions to the DE density 
respectively. In addition, the Hubble function H (normalized 
to 0.1 at the present time), the coupling function /3 C , the CDM 
particle mass m c (normalized to its present-day value), and the 
DE equation of state parameter w^, are included in the file. All 
the plots in Fig. [T] have been obtained using information which is 
included in these files. 

Other three files: 

motfeLCoDECS-growthf actor_CDM.dat 
moflfeLCoDECS-growthfactor _Baryons.dat 
mo<ie/_CoDECS-growthf actor _all.dat 

provide information on the evolution of the linear growth factor 
for CDM, baryons, and the total matter component, respectively, 
including the evolution of the quantity D+/a for these three 



different cases. Both plots in Fig. [2] were obtained using data 
included in these files. 

Finally, the linear matter power spectrum used to set up the initial 
conditions of the simulations and normalized to the amplitude at 
z = defined by erg for each specific model is also provided, both 
in terms of the dimensionless power spectrum A 2 (fc) and of the 
dimensional power spectrum P(k) = 2n 2 A 2 (k) / k 3 , in the files: 

morfeLCoDECS-linear _power.dat 

A3 Non-linear matter power spectra 

A series of links provides access to formatted files named: 

morfeLCoDECS-power_CDM_XXX.dat 
mtWe/_CoDECS-power_B aryons_XXX.dat 
m0deLCoDECS-power-all_XXX.dat 

for the nonlinear power spectra of CDM, baryons, and total matter, 
respectively, as extracted from the simulations at different redshifts. 
The number 'XXX' indicates the snapshot number, and a direct cor- 
respondence between snapshot number and redshift is provided for 
each link. As mentioned above, the power spectra have been com- 
puted at scales larger than the Nyquist frequency of each simulation 
by performing a CIC (Cloud-In-Cell) assignment to a grid of the 
same size of the PM grid, i.e. 1024 3 and 512 3 for the L-CoDECS and 
H-CoDECS sets of simulations, respectively. For Fourier modes be- 
yond the Nyquist frequency the power spectrum is computed with 
the folding method ( |Jenkins et al.|1998| |Colombi et al.|2009| > and 
then smoothly interpolated to the large-scale power in the region of 
overlap, in order to provide a single estimate of the nonlinear power 
up to modes where the shot-noise is found to be 10% of the mea- 
sured power, at which we cut our power spectrum estimate. Also in 
this case, the files include both the dimensionless and dimensional 
power spectra, A 2 (fe) and P(k), respectively. 

A4 Groups and SubGroups catalogs 

For all the CoDECS simulations, halos have been identified by 
means of a Friends-of-Friends (FoF) algorithm with linking length 
£ — 0.2 x d, where d is the mean interparticle separation. This 
procedure has been applied to the particle distribution in the sim- 
ulations for a large number of redshifts by linking the CDM parti- 
cles as primary tracers of the local mass density, and then attaching 
baryonic particles to the FoF group of their nearest neighbor. For 
each FoF halo we have computed a series of global properties as the 
total number of member particles, the total halo mass, the center- 
of-mass position and velocity. All this information, together with 
the number of gravitationally bound substructures identified within 
each single FoF halo, is included in the Groups files of the CoDECS 
database, that are grouped in archive files named: 

morfeLCoDECS-Groups JXXX.tar . 

Each archive contains 512 files for the L-CoDECS simulations 
and 128 files for the H-CoDECS simulations, corresponding to the 
individual files written by each processor in the different runs. 

The identification of the halo substructures is made by means 
of the SUBFIND algorithm ( Springel et al. 2001 1 and a large number 
of properties of the individual substructures identified in each FoF 
halo are computed. These include the location of the local gravi- 
tational potential minimum, the number of particles belonging to 



© 2011 RAS, MNRAS 000,[Tp8l 



The CoDECS project 17 



Simulation 


Box Size 


baryon particle 


CDM particle 


baryon particle 


CDM particle 


Gravitational 




suite 


L 


number 


number A^ c 


mass TYib 


mass tti c 


cnrtpnincT g 

AVJlLClllllJi tp 


Hydrodynamics 


L-CoDECS 


1 Gpc/h 


1024 3 
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Table Al. The numerical parameters of the two suites of runs of the CoDECS project: the collisionless large-scale suite L-CoDECS and the hydrodynamic 
small-scale suite H-CoDECS. 



the substructure, the total mass of the substructure, as well as its 
peculiar velocity and its specific angular momentum. The latter is 
defined as: 

5~] ■ rrii (ri x vA 
s = ^' A - (Al) 

where i runs over the particles bound to the subhalo, n is the posi- 
tion of the 2-th particle with respect to the potential minimum of the 
subhalo in physical coordinates, and Vi is its physical peculiar ve- 
locity. For the main substructure in each FoF halo, we also compute 
the virial radius -R200, defined as the radius of a sphere centered in 
the minimum of the gravitational potential of the substructure that 
encloses a mean density 200 times larger than the cosmic critical 
density, and the corresponding enclosed virial mass Mioo< These 
data are included in the SubGroups files directly available from 
the CoDECS web database, which are also grouped in archive files 
named: 

worfeLCoDECS-SubGroups JXXX.tar . 

Detailed information on the format and the units of the Groups 
and SubGroup data is provided in the CoDECS user's guide directly 
downloadable from the CoDECS web site. 
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